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Abstract. We present measurements of the triaxial dark matter halo shapes and alignment 
correlation functions in the Millennium and Millennium-2 dark matter A'^-body simulations. 
These two simulations allow us to measure the distributions of halo shapes down to 10% of 
the virial radius over a halo mass range of 6 X 10^ - 2 X 10^'' h-^MQ. We largely confirm 
previous results on the distributions of halo axis ratios as a function of halo mass, but we find 
that the median angle between halo major axes at different halo radii can vary by a factor 
of 2 between the Millennium-1 and 2 simulations because of the different mass resolution. 
Thus, error in the shape determinations from limited resolution is potentially degenerate 
with the misalignment of halo inner and outer shapes used to constrain Brightest Cluster 
Galaxy alignments in previous works. We also present simplifying parameterizations for the 
3-D halo-mass alignment correlation functions that are necessary ingredients for triaxial halo 
models of large-scale structure and models of galaxy intrinsic alignments as contaminants 
for cosmic shear surveys. We measure strong alignments between halos of all masses and 
the surrounding dark matter overdensities out to several tens of Mpc, in agreement 
with observed shear-galaxy and cluster shape correlations. We use these measurements to 
forecast the contribution to the weak lensing signal around galaxy clusters from correlated 
mass along the line-of-sight. For prolate clusters with major axes aligned with the line-of- 
sight the fraction of the weak lensing signal from mass external to the cluster can be twice 
that predicted if the excess halo alignment correlation is assumed to be zero. 
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1 Introduction 

The cold dark matter (CDM) model is now the paradigm for explaining the statistics of large- 
scale structure as well as galaxy formation and evolution. By using A^-body simulations it is 
possible to make detailed predictions of the cosmological distribution and evolution of CDM 
that can be compared with a wide array of observations. The large-scale statistics of CDM are 
a good fit to galaxy clustering and weak lensing observations while weak and strong lensing, 
galaxy velocity measurements, and galaxy cluster mass estimates are well-described by the 
predicted abundances and density profiles of virialized dark matter structures in A^-body 
simulations on scales larger than ~ 100 kpc. 

While the spherically averaged density profiles of dark matter halos first described 
by [1] (NFW) are good fits to both simulations and many observations, deviations from 
halo sphericity can be important for detailed comparisons of observations with the CDM 
model [2-4] . Measurements of halo mass density profiles for specific halos must account for 
asphericity to accurately compare with CDM models [5]. Halo triaxiality can also lead to 
biases in ensemble statistics such as the halo mass function and n-point correlations when used 
to constrain cosmological parameters [6] as well as measurements of the Hubble constant from 
clusters [7]. Accounting for cluster triaxiality is further important for detailed interpretation 
of lensing and Sunyaev-Zeldovich measurements [6]. 
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The details of the distributions of halo shapes and alignments also carry information 
about the hierarchichal formation of halos and the filamentary structure in which they are 
embedded [8-12]. In the current understanding of hierarchical structure formation the small 
mass halos form first from the anisotropic collapse of ellipsoidal overdensities in the mass dis- 
trubution. Larger halos grow both through fairly steady mass accretion of the surrounding 
dark matter and through "major mergers" with other halos. The "radial orbit instability" [13] 
can be important for determining the triaxial shapes of halos formed through steady mass ac- 
cretion, but it remains to be determined whether mergers with larger halos or the anisotropic 
accretion of mass from filaments play larger roles in the determination of halo shapes. Other 
recent work has shown that the triaxial density profiles of halos may be determined by the 
initial conditions of the mass-density perturbations [14]. 

Halos become more spherical with time as their constituent particles complete many 
orbits and undergo "violent relaxation" [15]. For a given halo mass and cosmological epoch, 
halos that formed more recently are also found to have smaller biases with respect to the dark 
matter distribution than halos that formed at earlier times [16, 17]. Correlations between halo 
shapes and formation times have also been found in A^-body simulations [18]. Together these 
results would seem to indicate that the halo formation time is a key parameter in determining 
the shapes and shape-dependent clustering of halos at given halo mass and epoch. However, 
the dependence of halo shape on formation time is nontrivial as many other halo properties, 
such as concentration, spin, and velocity and substructure statistics, have also been shown 
to be correlated with halo formation time and to infiuence the clustering of halos in ways 
that cannot be explained solely by the relationships between these parameters [18-22]. This 
complicated dependence of halo properties and clustering is referred to as "assembly bias" [20] 
and indicates that any model for the origins of halo shapes and orientation correlations should 
account for properties of the halo environment that have not yet been causally isolated. 

There are many recent attempts to analytically model triaxial halo shapes [23-27] based 
on the gravitational collapse of initally ellipsoidal mass overdensities. These models show 
good qualitative agreement with both observations and simulations, but simulations are still 
required to make precise predictions of the triaxial shape distributions and correlations in 
CDM because of the highly nonlinear gravitational growth over many decades in length scales. 

Predictions of large-scale clustering statistics using the halo model [e.g. 28] typically 
assume the spherically-averaged NFW profile but can be altered in detectable ways by halo 
triaxiality [29-31]. As shown by the pioneering work of [29], including triaxial halos in the 
halo model requires knowledge of the joint probability distributions of the halo axis lengths 
and orientations as functions of mass and redshift. Also, the "two-halo" term in the halo 
model requires a model for the correlations of halos as a function of the halo separation 
vector, and halo masses, axis lengths, and orientations. A major goal of the present paper is 
to provide these missing ingredients to the halo model. 

Also, the intrinsic galaxy alignments that contaminate weak lensing measurements are 
expected to be related to the intrinsic alignments of dark matter halo shapes. Previous groups 
have attempted to quantify the predicted galaxy intrinsic alignment signal by measuring 
projected halo shape correlations in A'^-body simulations and then comparing with observed 
projected shape correlations to calibrate the amplitude of the correlations [32, 33]. These 
groups further used the calibrated shape correlation function amplitudes as indications of 
the degree of misalignment between central galaxies and their parent dark matter halos. 
Because the distributions of projected triaxial shapes can be complicated to interpret [34], 
in this paper we focus on the measurements of 3D shapes and shape correlations in part 
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so that the modehng of ahgnments between galaxies and dark matter can be constrained 
robustly from observations. In this paper we compare measurements of triaxial halo shapes 
in the Millennium and Millennium-II simulations to both produce statistically significant 
measurements over a range of halo masses and to assess the convergence of the halo shape 
statistics with mass resolution. 

This paper is organized as follows. We briefly describe the A^-body simulations studied 
in section 2. Our definitions of halos and our methods for measuring shapes are described in 
section 3 and can be skipped by readers familiar with previous similar studies as well as those 
less interested in the details of our methods. In section 4 we update previously published halo 
shape distributions given our improved mass range and resolution and present new results 
on the orientations of halo shapes as a function of radius in section 4.2. Section 5 consists 
entirely of new measurements and analysis of the halo correlation functions binned in angles 
between the halo major axes and the major axis of one halo and the separation vector to 
another halo. We also present some projected halo alignment correlations that can be directly 
compared with previous work. We present simple parameterizations for the halo alignment 
correlations in section 6 and show how the measured alignment correlations can be applied 
to understand the bias in weak lensing cluster mass estimates from line-of-sight structures. 
We draw conclusions on the impact of halo alignments on several types of observations and 
describe future applications of our results in section 7. 

2 Simulations 

The Millennium^ simulation [35] solved for the positions and velocities of 10^*^ dark matter 
tracer particles, each of mass 8.6 X 10^ H-^Mq, in a cubic volume 500 h ^Mpc on a side. 
The Millennium-2^ simulation [36] used the same number of particles as the Millennium 
but in a volume 100 /i~^Mpc on a side so that the particle masses were 6.89 X 10^ h-^MQ. 
Both simulations were run the following cosmological parameters: 0^ = 0.25, JIa = 0.75, 
= 0.045, Us = 1, (Tg = 0.9, h = 0.73. This cosmology was chosen to be consistent with 
the WMAP3 [37] constraints. The improved cosmological constraints from WMAP7 [38] 
favor cTg w 0.8. While such a change in as will significantly alter the amplitudes of the halo 
correlation functions as well as the rate of nonlinear growth of structure, we do not expect 
large differences in the degree of halo triaxiality or alignments that cannot be understood 
with a simple time rescaling to match the linear growth functions as in [39]. 

3 Halo definitions 

In this section we describe our algorithm for identifying halos and determining their masses 
and triaxial shapes. 

We use the "Friends-of-Friends" (FoF) [40] catalogue generated for each simulation to 
initially identify the halos. For each halo, we read in all the particles in a sphere centered 
on the FoF center and with radius twice the FoF virial radius, r^ff , (as given in the FoF 
catalogue [36]). Next we redetermine the halo center by iteratively computing the center-of- 
mass of the particles within spheres whose radius is reduced by 2.5% [41] in each iteration. 
The iteration is stopped when no more than 500 particles remain in the sphere and the 
center-of-mass at this stage is taken as the new halo center. (We discard all FoF halos that 

^http : //www. mpa-garching .mpg.de/galf orm/millennium/ 
^http : //www . mpa-garching . mpg . de/galf orm/millennium- II/ 
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have fewer than 500 particles, which is consistent with the mass cuts we apply in Section 4.) 
We flag any halos whose new centers are more than 7% of the FoF virial radius away from 
the most-bound halo particle [42] as determined by the SubFind algorithm [43]. The offset 
between the most bound particle and the center of mass of a halo was previously used as a 
proxy for determining "unrelaxed" halos by [19, 42], where [19] also found that halos with 
large center offsets tend to have highly prolate shapes. By selecting only halos with center 
offsets less than 7%, we therefore intend to remove any bias in the shape distributions from 
unrelaxed halos. However, note that [42] used the center-of-mass of all particles within the 
virial radius while we determine the center-of-mass only from the 500 particles remaining 
at the end of the iteration for determining the halo center. Our flag for unrelaxed halos 
based on the center offsets is therefore liable to miss some halos that would be flagged by 
the algorithm in [42], but many of these are likely caught by our substructure flag described 
below. 

We then define updated halo masses and radii by determining the radius r of a sphere 
centered at the halo center where, 

^^^ = A..(.)peHt(.), (3.1) 

where M(< r) is the mass enclosed in the sphere of radius r and pcrit = 2>H'^{z) /SttG. We 
use two definitions for A^-: A2oo(-2) = 200 and [44-46] 

Ac(z) = ISvr^ + 82 - 1) - 39 {^ra{z) - if ■ (3.2) 

At the three simulation snapshots we will consider at z = 0, 0.5, and 1, the overdensity 
Ac(2;) = 94, 131, and 152. We define r/:^^{z) = r^\j:{z) using A^; = /^c{z) in equation 3.2 and 
use the notation r2oo to explicitly refer to radii computed at an overdensity of A200 (which 
is the overdensity often used to define the "virial radius"). We determine the radius that 
matches the target density by computing the enclosed density in expanding spheres; starting 
at and increasing the radii in increments of 0.005 r^^J^ until the density threshold is 

reached. 

We label the triaxial halo axis lengths as a(rsph) < ^(?^sph) < c(rsph). The axis lengths 
are defined to be functions of the spherical radius rgph from the halo center to allow for 
different halo shapes at different radii. For a given rgph, the axis lengths and orientations are 
defined to be the eigenvalues and eigenvectors, respectively, of the reduced inertia tensor of 
all the mass tracer particles contained within a sphere of radius r^ph centered on the halo, 

alo — part. 

/ii(rsph)= R^r^pl) ' ^'J = l'2,3, (3.3) 

where A'haio —part, is the number of particles in the halo, x^^i is the ith coordinate of the 
nth particle in the halo measured with respect to a fixed cartesian coordinate system. The 
normalization of the particle positions, Rn{r sph), is the elliptical radius of the nth particle 
defined in terms of the axis lengths, 

z'^ 

where x, y, z are defined in the principal-axis frame of the halo. Note that R is dimensionless 
because the axis lengths have physical units. The reduced inertia tensor gives less weight 
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to particles at the outer edges of the halo than the standard inertia tensor (without the 
l/i?^(rsph) weighting). We finally define the two 3D axis ratios as, 



a 

s = - 
c 



b 



(3.5) 



q = -. 

c 



Following previous analyses [9, 47, 48], we compute the reduced inertia tensor using an 
iterative algorithm that starts with all the particles within a sphere of radius rgph that is 
then deformed along the eigenvectors of the initial inertia tensor while keeping the volume 
within the (deformed) ellipsoid fixed. Because the axis lengths a(rsph), 6(rsph), and c(rsph) 
enter the definition of the inertia tensor explicity in equations 3.3 and 3.4, we rescaled the 
eigenvalues of the inertia tensor at each iteration by the quantity r^pYi/ {abc)^/'^ to impose 
the fixed volume constraint. Note that variations of the algorithm exist in the literature 
where the intermediate axis rather than the ellipsoidal volume is kept fixed between itera- 
tions. By keeping the enclosed volume fixed allows us to equate the shapes we measure at 
fixed spherical starting radius rgph with the shapes measured in [49] at fixed enclosed mean 
mass density. Particles that are outside of the ellipsoid boundary are dropped and any new 
particles within the ellipsoid are added in the subsequent recalculation of the reduced inertia 
tensor. The algorithm is terminated when either the axis ratios converge to a given toler- 
ance (max ((g' — qY /q^, {s' — s)^/s^) < 5 x 10~^) or the ellipsoid has fewer than 50 particles 
enclosed. This means the enclosed volume of the ellipsoid with axis ratios a, b, c is related 
to the spherical radius rgph from equation (3.4) as V = gT^pj^. Recently, [50] pointed out 
that our definition of the reduced inertia tensor may be more biased in recovering input tri- 
axial shapes than other possible inertia tensor definitions. However, the definition of inertia 
tensor must ultimately be chosen with the final application in mind. In this paper we aim 
to measure triaxial halo statistics that can be compared with other CDM halo predictions 
and also input into a halo model for galaxies. Because our reduced inertia tensor definition 
measures a shape of all enclosed particles, it is potentially more useful than, e.g., the shape 
of all particles in an ellipsoidal shell. The latter was used effectively by [49] to measure the 
triaxial density profile of halos, which we do not reconsider here. 

We restrict the inertia tensor calculation to those particles identified by the SubFind 
algorithm as gravitationally bound to the parent halo (denoted "SubFind halo-0" or SFO) [42] . 
That is, we consider only the halos that might host isolated galaxies or groups and clusters 
and discard all gravitationally bound satellite halos or other substructures. In addition, any 
particles that are neither gravitationally bound to the parent halo or in substructures are 
also discarded. We compute the reduced inertia tensor of SubFind halo-0 at 32 radii for each 
halo linearly spaced from 0.1 r^^ to rj\^. (Note these radii set the size of the sphere that is 
the starting point for the iterative inertia tensor calculation.) 

We discard all halos in the FoF catalogue whose centers, virial masses, or shapes do not 
converge in our iterative algorithms, but this is only 3% of halos in Millennium-1 and 1% 
of halos in Millennium-2 at z = (see appendix A for more detail on the halo counts with 
various cuts). We also discard all halos that have substructure mass fractions greater than 
10%, where the substructure mass is defined as the mass in all bound subhalos with masses 
greater than 1% of the virial mass as determined by the SubFind algorithm. The substructure 
mass cut further reduces the halo catalogues at z = by 24% and 15%, respectively, for the 
Millennium-1 and 2 simulations. The mass functions of the resulting halo catalogues are 
shown in figure 1. The vertical dotted and dashed lines in figure 1 indicate the mass bins 
that we consider in the remainder of the paper. 
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M200 (h-'M,„„) 



Figure 1. Comparison of the mass functions for all isolated ("SFO") halos in the simulations and 
those halos whose shapes were successfully measured ("cleaned"). The "cleaned" samples have cuts 
on the minimum number of particles in the halo (50), the shift in the halo centers and Virial masses 
between the FoF catalogue and our algorithms, the substructure mass fractions, and the sphericities 
(see text). The vertical dotted lines show the mass bins we use when describing the halo shape 
distributions while the vertical dashed lines show the mass bins for the halo correlation function 
measurements . 
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Our decisions to consider only the shapes of relaxed halos with substructures removed 
is motivated by our primary interests in the CDM predictions for halo shape statistics and 
in potential applications in modeling the shapes and orientations of central galaxies in the 
halo model. If we wanted to determine the elliptical halo shapes that could be detected in 
lensing measurements, we might instead measure the shapes of all FoF particles, perhaps 
in ellipsoidal shells rather than enclosed ellipsoids [50]. Although, our halo shape definition 
based on the reduced inertia tensor of all enclosed particles may in fact be more relevant for 
comparing with observable halo lensing signals due to the uncertainties of the line-of-sight 
projection [48]. In general, the choice of halo definitions and halo constituent particles will be 
application specific. We believe our choices are consistent with previous studies of CDM halo 
statistics and models for central galaxy alignments (which we explain further in Section 5.3). 
In Section 6.2 we use the shapes of SFO halos to estimate the bias in cluster weak lensing 
masses from correlated line-of-sight structures. Again, evaluating the validity of this choice 
will depend on the method used to select clusters and estimate their ellipticities. 

4 Halo shapes 

In this section we update the fitting functions for the joint distribution of the axis ratios first 
described in [49] to cover the larger mass and radius ranges available with the Millennium-2 
simulation [where we have particle masses of 6.89 x 10^ h~^MQ versus 6.2 x 10^ H^^Mq 49]. 
We also investigate the radial dependence of the axis ratios and the "twisting" of the halo 
axes with radius. We expect these new measurements to be particularly useful for modeling 
the halo mass distributions for lensing studies [e.g. 51] and for understanding the relative 
orientations of galaxies and their host halos. 

We do not update the triaxial density profile measurements of [49] because their high- 
resolution halo simulations already provided better resolution than in the large-volume sim- 
ulations we investigate here. However, we have confirmed that the triaxial density profile 
model of [49] is also an excellent fit to the halo shapes we measure. 

4.1 Axis ratio distributions 

The median axis ratios as functions of the elliptical halo radius are shown in figure 2. To 
avoid skewing the axis ratio distributions with halos that have poor shape determinations we 
have computed the medians after selecting only those halos with minor-to-major axis ratio 
s < 0.9 (which removes only 0.2% and 1% of halos in the Millennium-1 and 2 simulations, 
respectively, at z = 0). In all cases, the axis ratios increase with radius (meaning the halos 
become more spherical with increasing radius) and also increase with decreasing halo mass at 
fixed radius. The arrows in figure 2 show the "convergence radius," Tconvj defined so that the 
ratio, K = ireiax/icirc('"20o) = 7, where ireiax is the "collisional relaxation time" and tcirc('"20o) 
is the circular orbit timescale at r2oo [as defined in 41, 52]. Ref. [9] found that k = 7 marked 
the convergence radius where the axis ratios of individual relaxed halos agreed in different 
resolution runs of a quiescent halo in the Aquarius simulations of Milky Way-mass dark 
matter halos. For the intermediate mass bins where both simulations have good statistics, 
it is clear from figure 2 that the median axis ratios are still slightly misestimated for rconv 
defined at At = 7. This is likely due to the wide range of accretion histories of the halos 
in our catalogues that increase the scatter in the axis ratios at fixed halo radius. Because 
the Millennium-1 simulation has larger-mass particles, fewer substructures will be resolved 
and excised from the halos and will therefore contribute to the shape measurements. Also, 
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Figure 2. Median axis ratios of the triaxial lialo shapes as functions of elliptical radius from the 
center of the halo (normalized by r^i^ = r^^) in bins in log]^Q(Af2oo/^~^-^0)- The left panels show 
the ratio of the minor to major axis, s = ^, while the right panels show the ratio of the intermediate 
to major axis, q = K Note that all halos are restricted to have s < 0.9 to ensure reliable shape 
estimates. The axis ratios at redshifts 2 = 0, 0.5, 1 are shown from top to bottom. The arrows 
denote the "convergence radius" in each mass bin for each simulation where the ratio of the collisional 
relaxation time to the circular orbit timescale, k = 7. For a given simulation, lines arc omitted in 
mass bins that do not have complete sampling as indicated by the halo mass function. The vertical 
dotted lines denote the mean values of r2oo/r^^ for each redshift. 
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smaller-mass halos will tend to have larger mass-fractions of unresolved substructures and 
should therefore deviate more from the Millennium-2 results than the highest mass halos. 
This is exactly the trend we see in the discrepancies between the solid and dashed lines in 
figure 2. Because we exclude all halos with fewer than 50 particles within a spherical radius of 
O.lrAc, it is possible that our measurements of the axis ratio distributions are skewed at small 
radii in the lowest-mass bin for Millennium- 1 because we are selecting only high-concentration 
halos (the Millennium-2 mass bins all have significantly more than 50 particles). We have 
checked that increasing the minimum particle cut from 50 to 300 particles does not alter 
either the measurement of the convergence radius or the median axis ratios for radii above 
the convergence radius. Finally, we note that figure 3 in [48] showed that the definition of 
the inertia tensor has a noticeable effect on the distributions and median axis ratios, with 
the weighted inertia tensor yielding slightly larger axis ratios at all halo masses than the 
unweighted inertia tensor. However, the differences are smaller than the trends with mass 
and radius that we show in figure 2. 

Previous simulation studies have shown that more massive halos tend to have later 
formation times (as it takes longer to accrete mass and substructures for more massive 
halos) [16]. At a fixed time, massive halos then tend to be less relaxed than low-mass halos 
and are therefore less spherical and have lower concentrations [4, 53]. The decreases in the 
axis ratios with both increasing mass and redshift in figure 2 are consistent with this picture 
of the hierarchical formation of halos. 

The marginal distribution of the axis ratios at fixed mass is shown in figure 3. In figure 4 
we show the marginal distributions of the scaled axis ratio, 



in different mass bins, which are only weakly dependent on halo mass and redshift, where 
M^, is the mass where the r.m.s. density in top-hat spheres cj(M) is equal to the overdensity 
for spherical collapse, Ssc{z). The values of M*(z) for the three snapshots we consider are 
given in table 4.1. This scaling is motivated by [49] (their eq. 16), but we find we need a 



steeper mass dependence with index 0.0375[il(2:)]'^'^^ (= 0.030 at z = 0) rather than their 
index of 0.07[r2(2;)]^ '^ (= 0.027 at z = 0). While we study simulations with only one choice 
of cosmological parameters, [47] previously found that the dependence of the distributions 
of halo shapes on the amplitude of density perturbations, erg, was well described by the 
cosmology dependence of M* alone. 

We fit the marginal distributions for s shown in figure 4 with the same beta distribution 
(i.e. independent of mass). 





Table 1. Characteristic halo mass for our three snapshots. 



Redshift M*(/i-iM0) 



6.2 X 10 

0.51 1.3 X 10 

0.98 3.0 X 10 
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Figure 3. Marginal distributions of the axis ratios of the triaxial halo shapes at r^m and z = 0. s is 
the ratio of the minor to major axes and q is the ratio of the intermediate to major axes. The panels 
indicate bins in log]^g(Af2oo/^~^-^0)- 



where we find ~ 14.3 — 2.92; and fig ~ 10.4 — l.Sz as functions of redshift z. To describe 
the joint probability distribution for s and q we again follow [49] and plot the conditional 
distribution in figure 5, which then gives us p(s,q) = p{s)p{q\s). The conditional 

distributions for q naturally become narrower with increasing s because of the definition 
q > s. The distributions appear to be largely independent of halo mass as well as redshift 
(not shown). We again fit the conditional distributions of q given s with beta distributions 
as in equation 4.2. But, we first transform q onto the unit interval with the transformation 
q' = to have the same support as the beta distribution. We find s-dependent, but redshift 
independent, fit parameters for conditional distributions of q', aqi\g{s) « 72 — 230s + 222s^ 
and (3qi\s{s) « 1.6s~^-^. 

4.2 Halo morphologies and orientations 

The halo "triaxiality parameter" [47, 54], 

T={a^ -b'^)/{a'^ -c^), (4.3) 

is a convenient method of classifying the triaxial halo morphologies. We show the distribu- 
tions of the triaxiality parameter as a function of halo mass at z = in figure 6. A T < 0.33 
indicates an oblate ellipsoid, T > 0.66 indicates a prolate ellipsoid, while triaxial ellipsoids 
have 0.33 < T < 0.66. We confirm previous results that the halos become more prolate with 
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Figure 4. Marginal distributions of the scaled minor-to-major axis ratio s (defined in eq. 4.f) at r2oo 
and z = 0. The thick black line indicates the fit to the distribution of the axis ratios for all masses in 
the two simulations as described below equation 4.2. 



increasing halo mass [55] and also see that the shapes of the halos at small radii (0.25rA(.) are 
more prolate than the shapes at r a^. • The distribution of the triaxiality parameters at higher 
redshifts are qualitatively similar to figure 6, although the distributions at all masses and 
radii become slightly more skewed towards one (prolate) as redshift increases. We also note 
that for halos with masses > W^'^'^h~^ Mq it appears to be a good approximation to model 
the halos as prolate, so the shapes can be adequately described by the single minor-to-major 
axis ratio. We will make use of this approximation to simplify some measurements of the 
halo shape correlations in section 5. 

In figure 7 we plot the distributions of the angle between the halo major axis at O.lr^^ 
and the major axis at 32 logarithmically spaced larger radii out to r^^. The central lines of 
the boxes in figure 7 denote the median of the angles for all halos in each simulation while the 
box upper and lower edges denote the first and third quartiles of the angle distributions. It is 
interesting to see that the angle distributions for the two simulations are in good agreement 
in the largest mass bin (bottom panel) but are in larger disagreement as the halo mass 
decreases. We believe this is an artifact of the limited mass resolution in the Millennium-1 
simulation wherein some substructures are not fully resolved by the SubFind algorithm and 
are included as part of the main halo, leading to larger scatter in the halo shapes at different 
radii. In the higher resolution Millennium-2, more substructures are resolved and removed by 
SubFind, so that the parent halo has less twisting of the major axes as the radius increases. 
This is an important effect to note when modeling the alignments of BCG galaxies added to 
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Figure 5. Conditional distributions of intermediate-to-major axis ratios of the triaxial halo 
shapes given the minor-to-major axis ratio, s, at r2oo and z = 0. The panels indicate bins in 

logio(M2Oo//l"'M0). 

A^-body simulations as in [32, 56]. (We will return to this issue in section 5.3.) 

Because we have used the reduced inertia tensor (defined in equation 3.3) to compute 
the halo shapes, we expect that our measurements of the misalignment angles of halo major 
axes at different radii will be skewed more towards zero than if the un-reduced inertia tensor, 

-^halo — part. 

n=l 

were used to measure the shapes. This is because the reduced inertia tensor puts more weight 
on particles at small radii and is therefore more sensitive to the cumulative mass distribution 
in the halo. The choice of inertia tensor definition should therefore be kept in mind when 
comparing our results in figure 7 with other works. 

The wide distributions of angles in figure 7 is also remarkable. Because we have mapped 
all angles onto the interval (0, 90] degrees, the distributions of angles are highly skewed 
towards zero (indicated by the differences between means and medians). However, for radii 
> 0.5rA;, the distributions have large tails showing that ~25% of halos have outer major axes 
nearly perpendicular to the major axes at O.lrAc- To check that errors in the measurement 
of the halo orientation at ClrA^ have not skewed the distributions in figure 7, we recomputed 
the alignment angles relative to the axes at r^^ and confirmed that the relative alignment 
angles at different radii are robust to the reference radius. The wide distributions of alignment 
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Figure 6. Triaxiality parameter, T = (a^ — 6^)/(a^ — c^) at z = 0. The panels indicate bins 
in log]^o(M2oo/^^^A^0)- T < 0.33 indicates an oblate ellipsoid while T > 0.66 indicates a prolate 
ellipsoid. 



angles in figure 7 illustrate the necessity of simulations that can accurately resolve the inner 
regions of halos when modeling galaxy intrinsic alignments for analyzing weak lensing surveys. 
It is also possible that further investigations into these axis alignment distributions could 
show that the Milky Way is not unusual in having a distribution of satellites in a plane 
perpendicular to its disk [57]. 



5 Halo alignments 



To describe the alignments of the principal axes of different halos we measure the halo- 
halo and halo-mass 3D correlation functions binned in both the separation distance and the 
angles between either the halo major axes or the halo major axis and the separation vector 
to a tracer of the mass density (which we will call the "alignment angles"). Because the 
alignment effect is weak compared to the spatial clustering of the halos, we present most of 
our measurements in terms of the quantity, 
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Figure 7. Angle between the major axes at different radii within individual halos aX z — 1/2. Points 
indicate mean angles at a given radius while the boxes show the 25-75% central quantile range with 
the median of the distribution denoted by a horizontal line within each box. Each panel shows a 
different mass range with labels \ogn^{M20Qlh'~^ Mq). The halo radius on the abscissa is normalized 
by Tvir = (not r2oo)- 
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which gives the excess probabihty to find a halo or mass overdensity within a distance [r, r+dr] 
from a given halo and with the alignment angle within [9,9 + d9]. Note that, 

C(r) = / sm9d9^{r,9), (5.2) 
Jo 

so that the halo-halo or halo- mass correlation function ^(r) is the angle average of the align- 
ment correlation function ^(r, 0). When reporting our measurements we map all measured 
9 values onto the interval (0,7r/2) with the assumption that £,{r,9) = ^(r, vr — 9) due to the 
symmetry of the triaxial halo ellipsoids. We also ignore any potential dependence of the halo 
alignment correlation function on the orientations of the minor and intermediate axes. 

Another motivation for measuring the ratio defined in eq. (5.1) comes from [58] who 
recently showed that within the linear alignment model that assumes halo orientations align 
with the gradients of the large-scale gravitational potential, 

where 9p is the angle between the apparent major axis projected on the sky and the projected 
separation vector to the mass overdensity 6, and 7+ = (l — Sp) / (l + Sp) cos 29p with Sp the 
ratio of the projected minor and major axes. While this is defined in terms of the projected 
halo shape and orientation it gives us a reference for interpreting the radial and 9 dependence 
of S as defined in eq. (5.1) and is directly related to eq. (5.1) when the halo pairs are in the 
plane of the sky. 

We show measurements of the halo-halo and halo-mass position correlation functions 
in mass bins, ^(r;Mi,M2 = Mi), in figure 8 [see also 59]. The correlation functions are in 
good agreement between the two simulations in all overlapping mass bins from the minimum 
measured radii out to r ~ 20/i~^Mpc where the correlation functions in the Millennium- 
2 simulation begin to show artifacts due the imposition of the integral constraint in the 
100 /i~^Mpc simulation volume. In both rows of panels the separation of the correlation 
functions at different redshifts is clearly visible in the lowest mass bin due to the linear 
growth. For larger halo mass bins the increasing halo bias with redshift compensates for 
the decreasing growth and the correlations have similar amplitudes at all three measured 
redshifts. For the highest mass bin in the halo-halo correlation function (top-row, right 
panel), the effects of the nonlinear halo bias can be seen at small r where the amplitudes 
of the correlation functions are increasing with increasing redshift. At small-radii in the top 
row of panels in figure 8, the correlation functions turn towards zero due to halo exclusion 
(so the turn-over scale increases with increasing halo mass). In the bottom panels of figure 8 
the small-scale correlation instead show the correlations of the halos with the mass density 
interior to the halo (with the scale again set by the mass-dependent halo radii). We will use 
the correlation functions in figure 8 as references for assessing the effects of halo alignments 
in the next section. 

5.1 Angle-binned correlation functions 

In figure 9 we show measurements of C (r, cos 9; Mi, M2), with 9 the angle between the major 
axes of two halos in a pair and where the halo shape is measured at the r/\^. The error 
bars are estimated using fixed block bootstrap samples of the simulation volumes following 
the "marked-point bootstrap" method of [60]. In the marked-point bootstrap algorithm each 
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Figure 8. Halo-halo and halo-mass correlation functions in halo mass bins (miits are Mq/K) at 
redshifts 0, 0.5, and 1. The points where the halo-halo correlation functions turn over at small sepa- 
rations clearly the show the halo exclusion scales. As halo mass increases the large-scale correlations 
become more similar for all z because of the increasing halo bias with redshift. 



halo is assigned an array of "mark" values with each array entry equal to the number of halos 
(or mass tracers) in a given radius and angle bin relative to the first halo (so the mark array 
has length equal to the number of radius bins times the number of angle bins). We then 
divide each simulation volume into 64 equal-volume sub-cubes with fixed spatial locations 
and resample the sub-cubes 100 times with replacement to generate 100 realizations of the 
set of mark arrays in each simulation. For each set of resampled marks, we compute the 
correlation functions by summing the selected mark arrays for each halo, We then compute 
the ratio S(r, cos 6) for each resampling and determine errors by sorting the resulting values 
in each radius and angle bin to extract the 95% confidence intervals from the ordered arrays. 

In figure 9, we have plotted only the Millennium- 1 measurements in the three highest 
mass bins because the Millennium-2 measurements are noisy in these mass bins and there is 
no significant alignment signal in the two lower mass bins. However, there is a significant 
alignment signal that increases with halo mass with C~l.latr~l /i~^Mpc for the 
intermediate mass bins and the angle bin with < 9 < 22 degrees at z=0. The excess 
halo alignment correlations also increase with redshift giving an excess correlation of 1.25 at 
r ~ 1 /i~^Mpc and 9 < 22 degrees in the middle mass bin at z=l. In the CDM model we 
expect that the inner shapes of halos are in place at high redshift. As the evolution of the 
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Figure 9. Halo-halo C (r, 61; A/i, M2) versus r for bins in mass, \ogn^{Mi) and \ogn^{M2), and the 
angle, 9, between the major axes of two halo shapes at r = as measured in the Millennium- 
1 simulation at z = 0, 0.5, 1 from top to bottom. (The Millennium-2 measurements have larger 
uncertainties and are consistent with no alignments.) The shaded bands denote the 95% confidence 
intervals determined from 100 fixed block bootstrap samples. 
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halos proceeds, the filaments feeding mass onto the halos can re-form in new orientations 
causing shifts in the orientations of the halo shapes and decreasing the correlations with 
surrounding halo orientations. At high redshift there is also a greater range of length scales 
in the mass density perturbations collapsing at the same time, which also would support 
larger halo orientation correlations to larger halo separations. 

We show the halo-mass cross-correlation functions where 6 is now the angle between 
the major axis of a halo and the vector connecting the halo center to a mass tracer particle 
in figure 10. To compute the cross-correlations we selected a random subset of 10^ mass 
tracer particles in each simulation. There are very strong alignment signals for the halo-mass 
cross-correlation in all mass bins and both simulations are in good agreement. If there is any 
stochasticity between the alignment of halos and their surrounding mass distribution, then 
we would naturally expect the halo-mass excess alignment correlations to be much larger than 
the halo-halo excess alignment correlations. This is because in the latter case the random 
component of the halo alignments is included twice (once for each halo). Comparing the 
relative amplitudes in figure 9 and figure 10 would indicate a wide distribution in the angles 
between the halo major axes and the surrounding mass overdensities. 

It is also interesting to note in figure 10 that the alignment correlations are significant out 
to separations of several tens of megaparsecs for all measured halo masses, which is consistent 
with observations of galaxy intrinsic alignments [61, 62]. This is somewhat surprising however 
as the effects of the finite simulation box size are known to depress the position-position auto 
correlation functions on scales ~ 0.1 times the box size because the correlation functions 
must satisfy the integral constraint with the simulation volume. 

The excess correlations in figure 10 are also consistent with self-similar evolution at fixed 
M/M^{z). Using the M*(z) values from table 4.1 the lowest mass bin at z = 1, the second 
lowest mass bin at z = 0.5 and the middle mass bin at z = should represent roughly similar 
ratios of M/M^{z). Indeed we can see in figure 10 that the peaks in the excess correlations 
at, e.g., r ~ 1 /i~^Mpc have consistent amplitudes across the three redshift measurements. 

To better understand the peaks in the excess correlations around r = 1 /i~^Mpc in 
figure 10 we show similar correlation functions (for Millennium- 1 only) in figure 11 where the 
mass tracers are now SubFind-0 halos rather than mass particles as in figure 10. The separate 
rows of panels in figure 1 1 show mass bins in the halos whose alignments are measured while 
the columns of panels show mass bins in those halos that are used as mass tracers. Because 
we are now measuring the correlations of halo alignments with halo positions in figure 11 the 
correlation functions drop to zero at small radii just as in figure 9. For scales larger than 
the typical halo radii the excess correlations have similar shapes in figures 10 and 11 but the 
amplitudes are larger in figure 11, possibly because of the halo bias. At radii smaller than 
the peaks in the excess correlations in figure 11, the excess correlations quickly change sign 
(about 1) before hitting the halo exclusion radius. This is because the radius of the halo is 
larger along the halo major axis. So for radii < the typical halo major axis lengths there is 
greater excess probability to find a neighboring halo perpendicular to the major axis. The 
intermediate peaks around r = 1 /i~^Mpc in figure 10 are therefore where the two-halo excess 
correlations reach a maximum before changing sign. 

In appendix B we show the alignment correlation functions similar to those in figures 9, 
10, and 11 but with the halo shapes measured at ClrA^, rather than r^^. The amplitudes of 
the excess correlations are consistently smaller when the shapes are measured at O.lrA^ as 
expected if there is stochastic misalignment between the inner and outer halo shapes. 
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Figure 10. Halo-mass G {r,6; Mi) versus r for bins in mass, log]^Q(Afi), and the angle, 6, between 
the major axes of one halo (measured at r^^) and the position vector to a mass tracer. The shaded 
bands denote the 95% confidence intervals determined from 100 fixed block bootstrap samples. The 
dotted vertical lines show the maximum r/^^ value in each mass bin. 

- 19 - 



z = 




10"' 10" 10' 10- 10"' 10" 10' 10- 10-' 10" 10' 10= 10-' 10" 10' 10= 10-' 10" 10' 10= 

r (h"'Mpc) 




Figure 11. G {r,9; Mi, M2) where 9 is the angle between the major axis of one halo (measured at 
) and the separation vector to another halo. The rows of panels show the mass bins of the halos 
with the shape measurements. The columns of panels show the mass bins of the halos used as mass 
density tracers. The three plots show results for redshifts 2 = 0,0.5, 1 from top to bottom. 
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5.2 Axis ratio correlations 



Up to this point we liave considered only the correlations of halo orientations. We now 
consider how the halo-mass cross-correlation functions depend on the halo shapes as param- 
eterized by the minor-to-major axis ratio, s. The halo-mass correlations binned in both r 
and s, where s is determined from the halo shapes at r/\^, are shown in figure 12 normalized 
by the halo-mass correlation integrated over axis ratio s. The halo-mass correlations have a 
strong dependence on the minor-to- major axis ratio at all redshifts and in all but the most 
massive halo mass bin, where the statistics are weak at high redshift. At scales above the 
typical virial radius in each mass bin the ratio of '^(r, s)/^(r) is consistent with a constant 
value indicating a scale-independent excess halo bias that depends on the halo shape. Halos 
with large s values, corresponding to weak asphericity, are more biased than the total halo 
population while the most elliptical halos with small s values are less biased, as previously 
found by [18]. The excess halo shape bias in figure 12 appears to have only weak dependencies 
on halo mass and redshift. Also note that the bins with s G (0,0.5] break the general trend 
of the correlation function ratios with s in a few panels in figure 12 (z = middle mass bin 
and the lowest mass bin at all z). 

At scales smaller than the typical virial radii in each mass bin in figure 12, there is also 
a weak dependence on the axis ratio s, which is opposite that at large scales. [49] found 
that more prolate halos tended to have lower concentrations, which would boost the one-halo 
correlation function in the direction seen at small scales in figure 12. However, to limit the 
scope of the paper, we do not pursue this issue further. 

In principle, the halo shapes might also be correlated with their orientations. To test for 
such correlations we recompute the angle-binned halo- mass cross-correlations from section 5.1 
with each pair weighted by the quantity, 

1 - s2 

which is similar to the definition of the projected galaxy ellipticity given the projected axis 
ratio in weak lensing surveys. We compute the weighted correlation functions rather than 
computing a correlation function binned in both the axis ratios and orientation angles to 
increase the signal-to- noise ratio of the measurement. 

The resulting weighted correlation functions are shown in figure 13. We have first 
divided the weighted correlation function binned in alignment angle S^ui (r, 9) by the weighted 
correlation function integrated over alignment angles iwi^) to isolate the excess alignment 
correlation. We then divided by the identical ratio for the unweighted alignment correlations 
^(r, 9)/£^{r) (shown in figure 10). If the weight w was uncorrelated with the alignment angle 
6 then the "ratio of ratios" plotted in figure 13 would be equal to one for all r. The deviation 
from one indicates the degree to which halos with large values of w (i.e. highly aspherical halos 
with small s values) have excess alignment correlations relative to the entire halo population. 



We have plotted the measurements from the Millennium-2 simulation only in the two 
lowest mass bins in figure 13 to avoid cluttering the less-noisy Millennium-1 results in the 
three higher mass bins. For all mass bins the ratios are significantly different from one across 
the radius range from 0.1 to 10 /i^^Mpc, although the deviations are no more than 5%. Note 
that we truncated the range on the abscissa at 10 /i~^Mpc because the measurements are 
very noisy and consistent with no excess correlations on larger scales. The excess correlations 
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Figure 12. Ratio of tlie lialo-mass correlation functions binned in both separation, r, and minor-to- 
major axis ratio, s, to the correlation binned only in r at z = 0, 0.5, 1 from top to bottom. The halo 
shapes are measured at • The dotted vertical lines show the maximum r value in each mass bin. 
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Figure 13. Ratio of tlie lialo-mass alignment correlation functions weighted by 1 — s^/1 + to the 
alignment correlation with no weighting at z = 0, 0.5, 1 from top to bottom. The halo shapes are 
measm-ed at r/\^. 
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are also roughly constant over the plotted radius range and show very little change with halo 
mass or redshift. 

Both the tidal stretching of halos and anisotropic accretion of mass along filaments will 
tend to introduce both asphericity and halo alignments, so it is perhaps not surprising that 
the halo axis ratios should be correlated with their orientations. However, the mechanism 
causing the correlations of the halo alignments and shapes remains to be shown. 

5.3 Implications for measuring BCG alignments with their host halos 

Previous studies [most recently 32, 33, 56] have concentrated on the projected alignment cor- 
relation functions in A^-body simulations in order to more directly compare with observations. 
In this section we show that we can reproduce previous measurements of the projected align- 
ment correlation function by changing the definition of the inertia tensor used to compute 
halo orientations rather than invoking a stochastic misalignment of BCGs with their host 
halos as in the cited previous studies. This pertains to our earlier point that our reduced 
inertia tensor definition is appropriate for modeling the shapes and alignments of central 
galaxies. 

We compute Wp{r,p,9p) using an algorithm similar to that used to compute £,{r,6). To 
determine the projected angle Op between the major axis of the projected halo shape and the 
projected pair separation vector rp, we compute the 3D iterative inertia tensor as before to 
determine the particles belonging to a 3D ellipsoid but then compute the eigenvalue decom- 
position of the 2D inertia tensor of the final projected mass density. As in previous sections, 
the halo shapes are determined from SFO particles only (i.e. with bound sub-structures 
removed) and only for halos passing our convergence and relaxation criteria described in Sec- 
tion 3. From the projected alignment angles, we compute the anisotropic correlation function 
^(rp,n, 0p) and integrate the line-of-sight separation 11 over a slab of width 80 /i^^Mpc. 

The solid lines in figure 14 show the projected halo-mass cross-correlation Wp{rp,6p) 
binned in the angle Op between the 2D halo major axis measured at ta,, and the projected 
separation vector to the mass tracers. The dashed lines in figure 14 show the analogous 
measurement when our calculation of the inertia tensor in the iterative determination of the 
halo shape is modified from equation 3.3 to the unweighted inertia tensor lij oc ^XjXj, 
which is the inertia tensor used in [32] to define their "outer" halo shapes. Finally, the points 
in figure 14 show the measurements of Wp{rp,Op) using "red" galaxies in the SDSS by [32] 
(their figure 1 upper panels). To plot the points in our halo mass bins we have made a crude 
mapping of the magnitude bins used in [32] to halo mass according to the Halo Occupation 
Distribution model of [63] (their figure 3) for Luminous Red Galaxies in the SDSS catalogue. 
Note however that the "red" color selection in [32] may select non-central galaxies in the 
lower luminosity / halo mass bin. 

Figure 14 shows first that changing the inertia tensor from the unweighted to weighted 
definition causes the amplitude of the excess projected alignment correlation function to drop 
by a factor of ~ 2 in the largest mass bin (dashed to solid lines). Our improved match to the 
observations using the weighted inertia tensor is consistent with the result in [32] who found 
the best match to their observations using only the inner halo orientations (defined to enclose 
the equivalent BCG stellar mass) to model the BCGs because the weighted inertia tensor is 
also more sensitive to the inner halo shape than the unweighted inertia tensor. In contrast, 
[33, 56] defined halo shapes using the unweighted inertia tensor of all halo particles and then 
concluded a broad distribution of BCG alignment angles are necessary to fit their measured 
shape correlations. This suggests that (1) observations of Wp{rp,Op) (or similar statistics) 
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Figure 14. Projected alignment correlation function from MEl using two different algorithms for 
determining the outer SubFind-0 halo shapes. The "iterative" algorithm is described in Section 3 
and is used througout this paper. The "all SFO" method uses the unweighted inertia tensor of all 
SubFind-0 halo particles. The points show the measurements from [32] using "red" galaxies in the 
SDSS (the top panel in their figure 1.). The SDSS measurements were made in galaxy luminosity 
bins that are only approximately related to the halo mass bins plotted here. The 3 luminosity bins 
for the points from left to right are, -21 < Afo.i, < -20, -22 < Mq.^ < -21, -23 < Mo.u < -22. 



must be carefully compared with simulations (including a careful definition of "halo shapes" ) 
when drawing conclusions about halo alignments or alignments of BCGs with their parent 
halos, (2) the most useful definition of the "shape" of a halo is application specific and the 
reduced inertia tensor may be more useful for modeling the shapes of central galaxies (whereas 
the "shape" definitions advocated in [50] may be more relevant for the shapes derived from, 
e.g., weak lensing studies), (3) any misalignments between BCGs and the outer shapes of 
their parent halos might be sufficiently modeled by the misalignments of dark matter halo 
shapes at different radii (as shown in figure 7). Comparing the BCG alignment with the 
spatial distribution of satellite galaxies in groups [e.g. 64] can be a complimentary way to 
constrain BCG alignment angle distributions, which could potentially remove some of the 
ambiguity in modeling BCG orientations in simulations. 

6 Interpretation of halo alignment correlations 

In this section we consider how our measured halo alignment correlations could be incorpo- 
rated into the halo model [see, e.g., 28] for the dark matter two-point clustering statistics. 
A triaxial halo model was previously presented by [29], but lacked a well- motivated model 
for the correlations of halo shapes and orientations. Using a toy model for the triaxial halo 
correlations, [29] and [30] showed that including non-spherical halo shapes can lead to ob- 
servationally detectable changes in the predicted matter power spectrum with even more 
significant changes in the matter bispectrum. Perhaps even more critically for near-term 
observations, the halo model can be used to predict the contamination to galaxy weak lens- 
ing correlations from intrinsic alignments of galaxy shapes [65]. However, assuming galaxy 
shapes are correlated with the shapes of their parent halos, the correlations between triax- 
ial dark matter halo orientations and the surrounding mass overdensities (that would lens 
background galaxies) is a necessary input to any such halo model. 

Key components in the triaxial halo model are therefore the halo n-point correlations 
including dependencies on halo masses, ellipticities, and alignments. [29] call these the halo 
"seed" correlations, and we adopt their nomenclature. Using the shorthand from [29] for the 
joint probability distribution of Ni halos with positions Xj, masses Mj, axis ratios a,, and 
orientations Sj, 

p{l,...,Ni)= (6.1) 

p(xi, . . . ,XAr^,Mi, . . . ,MAr.,ai, . . . ,aAr.,ei, . . . ,£Ni), 

the two-point halo seed correlation function is defined by the relation, 

p(l, 2) =p(l)p(2) (1 + ^-^^(1,2)). (6.2) 

We now interpret the halo-halo alignment correlation functions from section 5.1 as the seed 
correlation functions of [29]. We may alternately choose to use the correlations measured at 
^Ac or ClrAc depending on the application. For example, the inner halo shapes are most 
relevant for modeling the alignments of central galaxies, while the outer halo shapes would 
be a better choice for modeling the alignments of clusters probed by the SZ effect or x-ray 
surveys [6]. 

The halo-mass correlations presented in section 5.1 can also be readily interpreted as 
seed correlations in the halo model with the assumption that the integral over the masses 
of the second halo in each pair has already been done. In figures 11 and 21 we present 
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measurements of the halo-mass ahgnment correlations when SubFind-0 halos are used to 
trace the mass density. These provide the seed correlations for modeling the halo-mass 
correlation in the halo model with the mass dependence of both halos in each pair explicit. 

6.1 Multipole decomposition of alignment correlations 

To make use of our measured alignment correlation functions in the halo model, it will be 
helpful to decompose the dependence on the alignment angle and the axis ratio s. 

We fit the angular dependence of the seed correlation functions with the separable 
model, 



where ^hx denotes the halo-mass or halo-halo correlation functions, Pn are the Legendre 
polynomials of degree n, and fn{f,s) parameterizes any radial and shape dependence that 
deviates from that in ^n{f)- Redshift dependence of the correlation functions and fit param- 
eters is implicit. Because f^^Pn{x)dx = the constraint that / £,{r, 6, s)d cos 9 = (,{r,s) is 
satisfied by equation 6.3. As pointed out by [29, 58], ^(r, 6*) also has the symmetry con- 
straints i^(r, 6) = ^(r, 6 + it) and ^(r, 9) = ^(r, —9), which are satisfied by our restriction of n 
to the positive even integers in equation 6.3. 

We have factored out ^h,x{f,x) in equation 6.3 to provide a simple fitting form for 
S{r,9) as presented in the previous section. By dividing ^f^^'^ by ^hx we explicitly remove 
the linear halo bias, leaving any residual stochastic or nonlinear bias to be absorbed in the 
fit parameters /n(r, s). 

Because Chsi^, s)/^h5{r) at large scales in figure 12 is roughly constant as a function of 
r at fixed s and halo mass and is also roughly constant as a function of halo mass at fixed r 
and s, we model. 



where bh is the usual linear halo bias as a function of halo mass, 6shape(S)-^) is a "shape 
bias" as a function of halo minor-to-major axis ratio and halo mass (see [21] for a similar bias 
factorization), and S^ss is the usual mass autocorrelation function. We estimate shape bias 
values by minimizing the mean squared error between the log of the halo-mass correlation 
functions binned in axis ratio s and the log of the halo-mass correlation integrated over 
s over the radius interval 3r^^ to 20 /i~^Mpc [16]. The estimated shape bias values are 
shown in figure 15 versus the peak-height threshold v{M,z) = ^^^^^^-^ where 5q,q{z) is the 
linear overdensity for spherical collapse [44] and ct(M, z) is the linear mass variance in top-hat 
spheres extrapolated to the redshift z. Plotting the shape bias as a function of ^'(M, z) allows 
us to combine the halo mass and redshift dependencies into a single variable and exposes 
the roughly linear increase in the shape bias for fixed s at all masses and redshifts. The 
coefficients of the linear regressions shown by the thick solid lines in figure 15 are given in 
table 2 where ishapels, J^) = hQ{s) + hi[s)u. Note that 6shape(s,-^) must obey the integral 
constraint (for given M and z), 




(6.3) 



"shape 



(6.4) 




(6.5) 



which must be imposed when using the values in table 2 in numerical models such as the 
halo model. 




Figure 15. Halo "shape bias" as defined in equation 6.4 as a function of the peak-height threshold 
i/(M, z) = (5sc(-z)/c(M, z) in bins in the minor-to-major axis ratio s. The colors show bins in the axis 
ratio while the point shapes denote the different rcdshifts where the halo-mass correlation functions 
were measured. The thick solid lines show linear regressions in each axis ratio bin while the shaded 
bands show standard errors on the regressions. We have excised the \ow-i/ values in the lowest bin 
in axis ratio (bottom red line and points) as measurements of these values are particularly noisy (see 
the top left panels in figure 12) 



Table 2. Linear regression coefficients for 6shapc(s, z^) 



Axis ratio bin 


^0 


bi 


(0, 0.5] 


0.21 


0.30 


(0.5, 0.6] 


0.33 


0.34 


(0.6, 0.7] 


0.68 


0.30 


(0.7, 0.9) 


1.13 


0.16 



- 28 - 



Integrating equation 6.3 over axis ratio, s, gives a model for the angle-binned alignment 
correlation functions presented in section 5, 



^seed 



1+ E fn{r)Pn{cOs{e)) 

n62Z+ 



where, 



fn{r) 



f Chxjr, x) fn{r, s) p{s) ds 
ihx{r) 



(6.6) 



(6.7) 



Fitting the excess correlations in figure 10 for /n(?') with n = 2,4 at each r yields small 
positive values of /4(r) that are consistent with zero at the 2.5-0" level for all halo mass 
bins and redshifts. Because we used only four bins in cos^ when measuring the alignment 
correlations we should expect to have a weak detection of higher-order angular dependencies. 
Below, we choose to focus only on the quadrupole term f2ir)- 

Finally, the s dependence of fnii", s) can be extracted using the weighted alignment 
correlation functions shown in figure 13. The weighted correlation functions with weights, 
w, are related to the binned correlation functions by, 



I Chx'^(^) ^) 'S) '"^(■s) p(w('S)) \dw/ds\ ds 



f wp{w) dw 



1+ Yl /-,nW^n(cOs(e)) 
rae2Z+ 



(6. 



Equations 6.7 and 6.8 can be rewritten with the aid of equation 6.4 to yield two integral 
equations for the shape-dependence of the functions /„ (r, s) , 



' fn{r, s) 6shape(s)p(5)ds = fn{r)^^^^ = ^^(r) 









1 u;(s) 
fnir,s) -r^6shapc(s)p(w^(s)) 
{w) 



dw 



ds 



ds = fw,n{r) 



(6.9) 
(6.10) 



The functions of r on the right-hand sides of equations 6.9 and 6.10 are the multipole mo- 
ments of i{r,0) and iwi^^O) normalized by the linear halo-mass correlation function. From 
equation 6.5 we see that equation 6.9 is satisified with /n(r, s) = F'^{r). However, if /n(r, s) 
is indeed independent of s then equation 6.10 implies, 



/ ^^6shape(s)p(^i'(s)) 
70 V^l 



dw 
ds 



ds 



Flir) 
Ff(r) 



constant 



Vn, r \ M. 



(6.11) 



We plot the ratio F2{r, M)/F"(r, M) for n = 2 measured from the halo- mass cross-correlation 
functions in figure 16. Note the two lowest mass bins show measurements in Millennium-2 
while the three higher mass bins show Millennium- 1 measurements. The ratio is roughly 
consistent with a constant value in each mass bin and redshift for scales larger than a few 
/i~^Mpc (note that measurements on scales larger than ~ 40/i~^Mpc are very noisy due to 
low statistics and are omitted from the plot). However, the ratios in figure 16 show large 
variations on scales between r2oo and ~ 3/i~^Mpc at z = 0, but become flatter as z increases 
in all but the least massive bin. It would appear consistent with our measurements then to 
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Figure 16. Ratios of weighted and unweighted quadrupole moments of the halo mass correlations. 
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neglect the s dependence of /2(?^, s) for r > 10/i^^Mpc at all z. At 2; = 1 and M > 10^^ h~^MQ 
the s depedendence of the quadrupole can be neglected for all r > r2oo- 

We note that the deviations of the ratios in figure 13 from unity does not necessarily 
require /„(r, s) to depend on halo shape, but could in principle be described by the shape 
dependence of the halo bias as described by the term bsha.pe{s, M, z). This can be seen 
explicitly by noting that the left-hand sides of equations 6.9 and 6.10 would both integrate to 
/„(r) if fn{r, s) = fn{f) and 6shape('S) were independent of s. But the correlation functions on 
the right-hand sides of these equations would also be equal in this case giving equal Legendre 
expansion coefficients fw,n{f) and fn{f) and therefore equal ratios in figure 13. However, 
because of the results in figure 16, we are led to conclude that both the shape dependence 
of the halo bias and the correlations of halo shapes and orientations could contribute to the 
behaviour in figure 13. 

It is also interesting to compare the fitted quadrupole coeficients to the quadrupole mo- 
ment of the linear mass power spectrum. These should be proportional under the assumption 
of the linear alignment (LA) model [58, 66-70] that the alignments of halo orientations fol- 
low the gradients of the large-scale (linear) gravitational potential. That is, because the LA 
model is based on the linear, Gaussian, mass density perturbations, the linear mass power 
spectrum is the only source for the halo-mass multipole moments. The amplitude of the 
LA model quadrupole is not determined by the theory, but deviations from the predicted 
scale-dependence should indicate where the LA model breaks down. It is expected that the 
LA model should be a poor approximation on scales of a few Mpc where it is known that 
filamentary structures dominant the shape of the gravitational potential around halos [8, 9]. 
In figure 17 we show the ratio of f2{''')ih&{f) / i2,&{f) ■, where. 



are the multipole moments of the linear mass power spectrum (which is largely insensitive 
the upper integration limit), and jn are the spherical Bessel functions of degree n. Note the 
lower integration limit in eq. 6.12 is set by the simulation box size, Lboxj which skews the 
large-scale shape of the multipoles compared to those in an infinite volume. Each panel of 
figure 17 is truncated at the lower range by r2oo corresponding to the most massive halo in 
each panel. The quadrupole moment ratios have steep negative slopes for all masses over the 
range r2oo < t < 40/i~^Mpc indicating significant differences from the LA model prediction 
(the Millennium-2 results in the two lowest mass panels are truncated at 1/4 of the box size, 
or 25/i"^Mpc). From this we can conclude that the alignments of halos with the surrounding 
mass overdensities over these length scales is sourced significantly by higher-order correlations 
in the nonlinear mass density field. 

6.2 Weak lensing cluster mass bias from correlated line-of-sight structures 

Galaxy cluster masses estimated from weak lensing are known to be contaminated by the 
lensing from line-of-sight structures external to the cluster because the lensing kernel is much 
larger than the size of a cluster [71, 72]. The line-of-sight lensing contamination is largest 
when the line-of-sight structures are spatially correlated with the cluster (so that there is on 
average more or less mass along the cluster line-of-sight than along a random line-of-sight). 
[72] showed that the contribution to the projected mass density in the vicinity of a cluster 
of mass M and orientation e from external line-of-sight structures within a slab of width L 




(6.12) 
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Figure 17. Estimates of the quadrupole moment, f2{r), of the halo-mass ahgnment correlation 
function, ^^^(r, 0) at r > r2oo- The estimates have been normalized by the ratio of the linear mass 
correlation function quadrupole and monopolc moments. If the halo alignments were sourced entirely 
by the quadrupole moment of the linear mass density field, then these ratios would be constant 
functions of r. The two lowest mass panels (left) show results from the Millennium-2 simulation while 
the three higher mass panels show the Millcnnium-1 results. 



centered on the cluster is given by, 



(Sext(r±|M,e)) = 2p^ / ^ dr. {1 + ^{r,e\M,e)) (6.13) 

J^rl^{M)-ri 

where ^ is the halo alignment cross-correlation function as shown in figure 10. An estimate 
for the total halo mass can be obtained by integrating the projected surface mass density over 
an aperture defined by W{r±). The contribution to the halo mass estimate from external 
mass along the line-of-sight is then, 

(Mext(M)) = J dr^_ W{r^) (Sext(r±|M, e)) . (6.14) 

Inserting eq. (6.6) into eq. (6.13), the contributions to the external mass are, 

(Mext(M, 7?)) = (M,„i(M)) + (Mcorr(M)) + (Malign(M, l9)) , (6.15) 

where, 

L 

{M^^,{M)) = 2pm j dv^W{v^) fl dry, (6.16) 
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and 

L 

{M^orr{M))=2pm f dv^W{v^) f] dr., CwW, (6.17) 

were considered by [72]. The new term we present here is, 

(Maiign(M,i9)) = 2pm [ dv ^ W{v^) [' dr,, Chsir) /2(r)P2(cos e), (6.18) 

and we have introduced the exphcit dependence on the angle between the halo major axis 
and the line-of-sight, •{}. Therefore, {McorriM)) is a measure of the weak lensing signal from 
correlated mass along the line-of-sight to a cluster that is isotropically distributed about the 
cluster. The new term (Maiign(M, ■!?)) on the other hand, is a measure of extra mass that is 
correlated with the cluster orientation. 

We compare (Mcorr(-^)) and (Maiign(M, ■!?)) with the total halo mass in figure 18 for 
two choices of M^(r_|_) presented in [72]. The left panel shows the fraction of the mass 
estimates from line-of-sight structures when integrating the projected surface mass density 
over a top-hat aperture with radius equal to r2oo- When ~ 0, (Maiign(M, i})) ?a (Mcorr(-^)) 
so that for triaxial clusters oriented along the line-of-sight the contamination from external 
line-of-sight mass is roughly twice as large as previous estimates. For ~ 7r/2, there is a 
decrement of mass along the line-of-sight so the total mass bias will be decreased relative to 
the isotropic case. In the right panel of figure 18 a compensated filter from [72] is used that 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

cos(t3) cos(#) 



Figure 18. Bias in the weak-lensing inferred halo mass from correlated structures along the line 
of sight as a function of the cosine of the angle of the halo major axis with respect to the line-of- 
sight. The left panel assumes a top-hat filter is used to integrate the projected surface mass density to 
estimate the total halo mass. The right panel assumes a compensated filter that is chosen to maximize 
the signal-to-noise ratio of the halo mass estimate for a spherically symmetric NFW halo. 

is chosen to maximize the signal-to-noise ratio of the halo mass estimate when the halo has a 
spherically symmetric NFW density profile. Both (Mcorr(-^)) and (Maiign(M, -3)) are reduced 
to negligible levels for current cluster weak lensing studies. However, the compensated filter 
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requires finding a good guess for the true halo mass and concentration. Any errors in this 
guess win result in increased mass bias terms. 

As we mentioned in Section 4, the predicted bias in cluster weak lensing mass estimates 
could change with a different choice of inertia tensor and halo constituent particles used to 
measure the halo-mass alignment correlation function. These choices will depend on the clus- 
ter selection method and the estimator for cluster ellipticity when applying our predictions 
to interpret observations. For example, our results as presented may be useful if the BCG 
shape and orientation is used as a proxy for the cluster orientation while a different shape 
estimator may be more useful to compare with x-ray or Sunyaev-Zel'dovich morphologies. 
We leave the determination of shape estimators for these applications to future work. 

7 Conclusions 

We have measured the triaxial shapes and the alignment correlations of dark matter halos 
over wide ranges in halo radius (O.lrAc - '''Ac) ^^'^ halo mass (10^*^ - 2 x 10^"^ Mq) by 
combining measurements in two similuations with different volumes and mass resolutions. 
While we largely find good agreement between the two simulations in the mass ranges where 
they overlap, there are some important discrepancies. First, near the lower mass limit for 
the Millennium simulation (with particle mass 8.6 x 10^ h~^MQ) the distributions of axis 
ratios are slightly skewed towards more spherical halos than in the Millennium-2 simulation 
(with particle mass 6.885 x 10^ h~^MQ). Second, the distributions of angles between the halo 
major axes at small radii and large radii are narrower and closer to zero in the Millennium-2 
simulation indicating less misalignment between halo shapes at different radii. 

We draw several conclusions based on the distributions of halo shapes and the agreement 
or discrepancies between our two simulations: 

• The mass resolution of the A^-body simulation and the method for computing the in- 
ertia tensor of halos can significantly alter the alignments of halo axes as a function of 
halo radius as well as the correlation of halo orientations with the surrounding mass 
distribution. Because we are using halo particles that are identified to be graviation- 
ally bound but not part of any bound substructures, we have speculated that more 
substructures are resolved and removed in the higher-resolution simulation. In general, 
we would expect unresolved substructures to both make the inferred halo shapes more 
spherical and to randomize the orientations of the halo shapes at different radii, in 
agreement with our results. Constraints on the alignments of BCGs with their host 
halos should account for these issues. 

• The mean angles between the halo major axes at small and large radii are typically 
~ 20°, which is significantly less than other claims in the literature [33, 73]. However, 
there is very large scatter in the distributions of alignment angles so that ~ 25% of 
halos of all measured masses have nearly perpendicular major axes at small and large 
radii. This is a somewhat surprising result, indicating that even if we assume central 
galaxies trace the inner shape of their parent halos, we should frequently expect large 
misalignments with tracers of the outer halo shape such as the satellite distribution or 
gravitational lensing measurements. The large misalignments of inner and outer halo 
shapes therefore confound constraints on the dynamical accretion of satellite galaxies 
based on their spatial distributions. 
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• Halo shapes become less spherical with increasing halo mass and redshift and decreasing 
spherical halo radius, confirming the results of previous simulation studies[47, 49], but 
extending to larger mass and radius ranges. However, simulations including baryons [74] 
indicate the axis ratios can increase by as much as 0.4 at small fractions of the virial 
radius (O.lrigo), with systematic increases in the axis ratio for all radii. This remains 
an important issue to consider when confronting our measured halo shape distributions 
with observations. As shown in figure 6, halos of all masses become systematically 
more prolate for decreasing spherical halo radius (equivalent to constant enclosed mass 
density in our measurements). For halos with virial masses < W^^h-^MQ the shapes 
at the virial radius tend to be oblate while remaining strongly prolate in the interior 
of the halo. For larger halo masses the halo shapes are predominantly prolate at all 
radii, indicating that for many applications high-mass halo shapes can be reasonably 
modeled with only the value of the minor-to-major axis ratio. 

To study the correlations in both the shapes and orientations of halos with the large-scale 
structure, we computed several types of correlation functions binning in the minor-to-major 
axis ratios of the halos or the angles between the major axes of two halos or the major axis 
of one halo and the separation vector to a mass density tracer. The conclusions we draw 
directly from these measured correlations include: 

• The excess halo alignment correlations are strongly increasing functions of halo mass, in 
agreement with the paradigm that high-mass halos are both younger than lower-mass 
halos and more biased (i.e. are located at the connecting points between filaments). 

• The excess halo-mass alignment correlations are significant to several tens of mega- 
parsecs, again in agreement with measurements of shear-mass correlations [62] and 
group and cluster alignment correlations [75, 76]. 

• The excess alignment correlations are both much larger and more significant in the 
halo-mass correlation functions (figure 10) than in the halo-halo correlation functions 
(figure 9). We speculate that there is a large stochastic component to the alignments 
of halo major axes with the surrounding mass distribution, which serves to further 
decrease the amplitude of the halo-halo alignment correlations because the alignments 
of both halos in each correlated pair must be averaged over. 

• Halos with small minor-to-major axis ratios are less biased (i.e. have smaller correlation 
function amplitudes) than the average halo population at fixed mass and redshift. 
This is consistent with previously established relations between the age and bias of 
halos [16, 20, 21, 77] if we understand that more spherical halos typically formed earlier 
(at fixed halo mass and redshift). 

• Our measured halo-mass alignment correlation function can be adequately described by 
only a quadrupolar angular dependence, although we do detect a statistically significant 
nonzero coefficient of the fourth order multipole for higher mass halos. We are limited 
in this measurement by our use of only four angular bins in measuring the alignment 
correlations. 

• The alignment quadrupole has a shape that is significantly different from the quadrupole 
of the linear matter correlation function (see figure 17), which is the predicted shape 
in the linear alignment model [58]. 
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As one direct application of our measured halo-mass alignment correlation functions, 
we updated forecasts for the contamination to galaxy cluster mass estimates from weak 
gravitational lensing observations due to correlated (but not gravitationally bound) structures 
along the line-of-sight to the cluster. Previous studies such as [72] showed that the halo-mass 
correlation function implies average lensing contamination for clusters that can be ~ 10 — 20% 
of the cluster mass when using a top-hat aperture to estimate the mass. We introduced a 
new contamination term that is correlated with the orientation of the prolate cluster with 
respect to the line-of-sight and that can give additional mass biases similar in magnitude to 
those estimated in [72]. It has already been shown that the bias in cluster mass estimates 
from line-of-sight contamination is significant for constraining cosmological parameters from 
cluster counts [78] and our new bias term will only increase the error in inferred parameters. 
However, if the orientation of the cluster with respect to the line-of-sight can be independently 
observationally determined (e.g. [79]), then it may be possible to reduce the scatter in weak 
lensing mass estimates by subtracting the mean bias from correlated structures as predicted 
by the halo-mass alignment correlations we have measured. It is also possible, but remains 
to be investigated, whether the mass correlated with cluster orientations is a significant bias 
for strong lensing measurements in clusters (see e.g. [80]). 

Our alignment correlation function measurements will also be useful for building halo 
models for intrinsic alignments of galaxies as contamination in cosmic shear surveys [65, 81]. 
The so-called "intrinsic-intrinsic" (II) contamination for cosmic shear is best modeled using 
our halo-halo alignment correlations based on the inner halo shapes as shown in figure 19. 
Note that such an II model would be built by first populating central galaxies with 3D 
shapes that match the 3D shapes of their parent halos and then performing the line-of-sight 
projection to predict the 2D cosmic shear contamination. This would further reduce the 
amplitude of the alignment signal, so we expect any predicted II models to be quite small 
on scales larger than a few megaparsecs. The so-called "galaxy-intrinsic" (GI) cosmic shear 
contamination term is likewise best modeled with the halo-mass alignment correlations shown 
in figure 20 or 21 and modeled in equation 6.6. 

Because our measured alignment correlations have halo mass dependence that is roughly 
self-similar (see the discussion in section 5.1), and because we do not expect the scaled 
distributions of halo shapes to be strong functions of cosmology [47] , it may become possible 
to include the cosmology dependence in models of galaxy intrinsic alignments simply through 
the dependence of the alignment correlations on the characteristic mass scale M^:(z) and the 
matter power spectrum. However, one should keep in mind that halo occupation statistics and 
galaxy properties influencing their observational selection probably depend on a number of 
halo properties that are related by the as yet un-modeled "assembly bias" [16, 18, 20, 21, 82], 
which will confuse any cosmological dependencies in the intrinsic alignments in cosimc shear 
surveys. 

Because the halo shape statistics at small halo radii (where the galaxy alignments are 
often modeled) are slow to converge with mass resolution, galaxy intrinsic alignment models 
may be more accurate when using a halo model calibrated from high-resolution simulations 
than when populating halos with galaxies in a large- volume simulation (which will generally 
have lower mass resolution due to computational constraints). We will pursue this approach 
in a forthcoming publication. 

Finally, the alignment between galaxies and dark matter is an additional complication 
when comparing our results with observations [83, 84]. Both the effects of gas on the shapes 
of the dark matter halos and the relative orientations of galaxy shapes and their parent halos 
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remain promising avenues for future research. 
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A Halo counts 

We present counts of the numbers of halos with shape measurements in both our simulations 
in tables 3 (for Millennium-1) and 4 (for Millennium-2). The three columns in each table 
show counts for the three simulation snapshots considered in this paper. The rows in the 
tables show the numbers of halos when different cuts are applied to the halo shape catalogues. 



- 41 - 



Table 3. Halo counts for Millennium- 1 



Description 


z = U 


z = 0.5 


z = 1 


No cuts 


20,350,503 


20,441,331 


20,444,688 


In shape mass bins 


16,217,712 


17,078,569 


16,507,354 


Quality cuts in shape mass bins 


15,622,670 


15,930,109 


14,815,635 


Substructure cuts in shape mass bins 


12,372,629 


12,124,709 


10,979,741 


Axis ratio cuts in shape mass bins 


12,349,801 


12,115,975 


10,976,309 


In corr. mass bins 


15,014,590 


15,846,000 


15,173,960 


All cuts in corr. mass bins 


11,439,534 


11,285,980 


10,148,373 



Table 4. Halo counts for Millennium-2 



Description 


z = 


z = 0.5 


z = l 


No cuts 


11,073,314 


15,921,384 


12,200,000 


In shape mass bins 


8,386,084 


9,430,983 


9,559,984 


Quality cuts in shape mass bins 


8,307,617 


9,254,432 


9,266,159 


Substructure cuts in shape mass bins 


7,147,623 


7,626,275 


7,313,926 


Axis ratio cuts in shape mass bins 


7,075,515 


7,588,251 


7,295,991 


In corr. mass bins 


3,692,048 


4,156,360 


4,146,161 


All cuts in corr. mass bins 


3,067,434 


3,287,135 


3,103,102 



The first row, labeled "No cuts," gives the total numbers of halos with shape measurements 
irrespective of any cuts, except that the Millennium-2 halos with virial masses less than 
9.3 X 10^ h~^MQ were not considered. The second rows in each table, labeled "In shape mass 
bins," give the halo counts in all the halo mass bins depicted by the vertical dotted lines in 
figure 1. The third rows, labeled "Quality cuts in shape mass bins," are the remaining counts 
in the mass bins after halos are discarded that had unconverged centers or virial masses. The 
fourth rows, labeled "Substructure cuts in shape mass bins," give the halo counts when an 
halos with more than 10% mass in substructures are discarded, in addition to the previously 
applied cuts. The fifth rows, labeled "Axis ratio cuts in shape mass bins," give the remaining 
halo counts when an additional cut is applied that the minor-to-major axis ratio s < 0.9. The 
sixth rows, labeled "In corr. mass bins," give the counts in the mass bins used for measuring 
the halo correlation functions as depicted by the vertical dashed lines in figure 1. The final 
seventh rows, labeled "All cuts in corr. mass bins," give the counts in the correlation function 
mass bins when all the previously listed cuts are applied. 

B Halo alignment correlation supplementary plots 

In figure 19 we show the halo shape-shape correlations similar to figure 9 but with shapes 
measured at O.lrA^,. There is no detectable alignment of inner halo shapes at z = 0, and only 
a moderate alignment signal at z = 1. This indicates that the misalignment of the inner and 
outer halo shapes is not strongly correlated with the intra-halo alignments. 

The halo-mass excess alignment correlation functions using the halo shapes at 0.1r/\^ 
are shown in figure 20. Contrary to the halo-halo alignment correlations with the inner halo 
shapes, there is still a large alignment signal for the halo-mass correlations. The peak ampli- 



- 42 - 



z = 



1.50 - 

1,25 - 

1.00 - 

0,75 - 

_ 0.50 - 

+ 1.25 - 

S ' "0 " 
CD 

^0.75 - 

+ 

1.50 - 
1,25 

1.00 - 

0,75 - 

0.50 - 







M; (11.8, 


2.1) 




1 

H 


M; (12.7, 


3,5) 






M; (13.5. 1 


4,3) 



























































































1 
































F 


ai..^ 


1: (12,7 
























S 












































vm 




1 


















1 




^ 1 




















f 


























1 



Angle bins [deg.] 

(0,22,5] 
(22,5.451 
(45,67,5] 
(67,5,90) 



10"' 10° 10' 10- 10"' 



10" 10' 
r {l)''Mpc) 



10- 10-' 10° 10' 10- 



z = 0.5 



1,50 - 

1,25 - 

1.00 - 

0,75 - 

_ 0.50 - 

-I- 1.25 - 

S - 

CD 

vSo.75 - 

+ 

— 0,50 - 

1.50 - 

1,25 

1.00 - 

0,75 - 

0,50 - 



r {l)''Mpc) 
Z= 1 



1,50 - 

1,25 - 

1,00 - 

0,75 - 

_ 0,50 
^ 1,50 
m!? 

-I- 1.25 - 

S - 

CD 

v£o.75 - 

+ 

— 0,50 - 
1.50 - 

1,25 

1.00 - 

0,75 - 

0,50 - 



10-' 10° 10' 10- 10-' 10" 10' 10- 10-' 10° 10' 10° 







M; (11.8. 


2.7) 




M: (12.7, 13.5) 






M: (13.5. 14.3) 
















































1 






=1 



























































































r 


















































1 

2 




































































HI 



(22,5.45] 
(45,67,5] 
(67,5,90) 





L 


M: (11.8, 


2,7) 




1 


M: (12.7, 


3.5) 




M: (13.5. 


4,3) 






i 


















i 




J 








=i 
















r 






































fcl 




1 
























s 

1 


















m 






















f 




-\ 


















1 


|[ 




1 — 1 





































































Angle bins [deg.] 

y (0.22.5] 
(22,5,45] 
(45.67,5] 
(67,5.90) 



10-' 10° 10' 10° 10-' 10" 10' 10° 10-' 



10° 10' 



r (])~'Mpc) 



Figure 19. 6 (r, 9; Mi, M2) as in figure 9 except with shapes measured at r — O-lrA^ 
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tudes of the excess alignments are reduced by about 20% relative to the halo-mass alignments 
shown in figure 10 using the halo shapes at r^^. The shapes of the excess correlations are also 
changed relative to figure 10, with a reduced "bump" on scales just larger than the typical 
virial radii in each mass bin. 

Finally, in figure 21 we show the halo- mass alignment correlations where SubFind-0 
halos are used as mass density tracers and the shapes of halos are measured at O.lrA^.. 
Again, the amplitudes of the peak alignment correlations are reduced by ~ 20% relative to 
those in figure 11 based on the halo shapes at r/\^. We expect the alignment correlations 
shown in figures 20 and 21 to be most useful for modeling the alignments of central galaxy 
shapes, which should more closely trace the inner rather than the outer shapes of their parent 
halos. 
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Figure 20. G {r,0; Mi) with 9 the angle between ttie major axes of one lialo (measured at CItaJ 
and tlie position vector to a mass tracer. Tire shaded bands denote the 95% confidence intervals 
determined from 100 fixed block bootstrap samples. 
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Figure 21. G {r,9; Mi, M2) where 9 is the angle between the major axis of one halo (measured at 
O.lrAj and the separation vector to another halo. The rows of panels show the mass bins of the halos 
with the shape measurements. The columns of panels show the mass bins of the halos used as mass 
density tracers. The three plots show results for redshifts 2 = 0,0.5, 1 from top to bottom. 
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